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Abstract 

A theory is presented for the modification of bandgaps in atomically thin boron nitride (BN) by attractive interactions 
mediated through phonons in a polarizable substrate, or in the BN plane. Gap equations are solved, and gap 
enhancements are found to range up to 70% for dimensionless electron-phonon coupling A. = 1 , indicating that a 
proportion of the measured BN bandgap may have a phonon origin. 
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Background 

The need for bandgaps in graphene on electronvolt scales 
has led to a number of proposals, such as the use of bilayer 
graphene [1], creation of nanoribbons [2], and manipu- 
lation through substrates [3,4]. Recently, it has become 
possible to manipulate atomically thin layers of boron 
nitride (BN) and other materials with structure similar to 
graphene [5]. This may lead to a complimentary method 
of manipulating bandgaps to make digital transistors. 

In low dimensional materials, strong effective electron- 
electron interactions can be induced via an interaction 
between electrons confined to a plane and phonons in a 
polarizable neighboring layer [6]. The theory has shown 
that similar interactions account for the transport prop- 
erties of graphene on polarizable substrates [7] and that 
sandwiching graphene between polarisable superstrates 
and gap opening substrates can cause gap enhancement 
[8] . This paper examines similar gap changes in atomically 
thin BN due to interactions mediated through substrates. 

Methods 

Atomically thick hexagonal BN (h-BN) has similar chem- 
istry to graphene: bonding occurs through sp2 hybridiza- 
tion, and electrons with energies close to the chemical 
potential are in unhybridized tt orbitals [9]. A key dif- 
ference is that the electronic charge is not completely 
screened by the sp2 hybridization, shifting tt orbitals by 
An = +A on N sites, and — A on ^ sites. This shift is the 
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dominant cause of a gap of order 2 A. Tight-binding fits to 
results from ab-initio simulations of monolayer BN have 
established the hopping t=233 eV [10], with an estimate of 
A = 1.96 eV=0.84 1, The experiments indicate larger gaps: 
bulk h-BN has 5.971 eV [11], and monolayer h-BN has a 
gap of 5.56 eV [12] corresponding to A = 2.78 eV=1.20 t 
There is significant variation in phonon energies, hQq, 
in h-BN [13]. Longitudinal acoustic (LA) phonon ener- 
gies range up to around 140 meV at the M point, and 
transverse acoustic (TA) modes to around 110 meV at 
the K point. Optical phonon energies range between 160 
and 200 meV. Coupling, fn(ift), between electrons and 
phonons in either a polarisable substrate or the BN mono- 
layer, is possible and the corresponding Hamiltonian is as 
follows: 



H 



t {a 

{n,n')a ni 
hQ(Nm + 1/2) + Y ^nnn 



(1) 
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The Hamiltonian terms are shown schematically in 
Figure 1 (a), ali^j creates electrons of spin a on B sites and 
c^,^ on N sites. Vectors n are to atoms in the monolayer, 
and m to atoms in the substrate. Nm and yina are the num- 
ber operators for phonons and electrons respectively. §^ is 
the atomic displacement. The Hamiltonian is also approx- 
imately valid for interactions in the plane, and Figure 1 (b) 
shows the forces on ions from an increase in electron den- 
sity at a B site. The largest forces are on the near-neighbor 
sites, so that the effective interaction is mainly site diag- 
onal (electrons on A sites self-interact through phonons 
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Figure 1 Bn substrate system and interactions in a monolayer of 

BN. (a) BN substrate system annotated with interactions. 
Electron-plionon interactions between tine BN layer and substrate are 
poorly screened, and large interactions of strength fn(tn) are possible. 
Ions in the substrate oscillate with frequency Q. N sites have energy 
H-A and B sites - A, opening a gap. The attractive phonon-mediated 
electronic interaction / binds electrons onto the same site, effectively 
enhancing the gap. (b) Interactions in a monolayer of BN. Red circles 
represent N atoms and black circles B atoms. Light blue arrows 
represent distortions expected from an excess of charge on the site 
labelled e. 



modulation of the electron-phonon interaction due to 
incommensurability of the substrate, which was estimated 
at around ±8% of the average value [8] . 

Results and discussion 

The low order perturbation theory is applicable for low 
phonon frequency and weak coupling. A set of gap 
equations was derived by symmetrizing the self energy, 

^ ihconil - Z^) + Ai 0 



The local approximation used here is a good starting 
point because the modulated potential A is large, and 
electrons are well localized. Off-diagonal terms do not 
feature in the lowest order perturbation theory for the 
Holstein model since the interaction is site diagonal. 
is the quasi-particle weight and An is the gap function. 
For bosonic quantities, hcos = 27tkBTs and for fermions, 
hcoyi = 27TkBT{n + 1/2). T is the temperature and n and s 
are integers. 

The full Green function can be established using 
Dysons equation G~^{k,icL>n) = GQ^ikJcOn) — ^{icon), 
leading to. 



Z^ihcon - A - 

0Ar 



Z^ihcDn + A + Af 



(2) 



Substituting the expression for the Green function into 
the lowest order contribution to the self energy. 



Vbz 



Gii (k - q, i(^n-s) 



X (d^^{q,cOs)-2d^^{0,0)). (3) 

Here, the phonon propagator, (Iq'^ (q, cog) = SijQp'/{Qp' + 
co^y, hence, there are no off-diagonal elements of the low- 
est order self energy. The use of a single averaged Q and k 
here is consistent with the mean-field approximation. At 
half-filling, it is reasonable to assume that Xa = Xbso that 
A"^ = and = Z^. This leads to the gap equations. 



-tXkBTj2 j 



Die) A'„_^(do(icos)-2) 



(4) 



on B sites and vice versa). The diagram indicates that 
the strongest interaction is between electrons and optical 
phonon modes. 

For simplicity, the Holstein electron-phonon interac- 
tion was used,^(m) oc S^^n'^ which qualitatively captures 
the physics. There may be quantitative changes to the 
results for longer range Frohlich interactions and from 



^"-^ a>„ ^l'''n^col_,zl_, + Ats + ^'' 

(5) 

where the full gap is A^ = A^+ A. The density of states for 
a tight binding hexagonal lattice in the absence of a gap, 
D{e)y has the form given in reference [14]. The equations 
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may be solved self-consistently by performing a truncated 
sum on Matsubara frequencies. 

Gap and quasi-particle weight functions only have a 
weak Matsubara frequency dependence (< 0.3% for k = 
1, k^T = hQ = 0.01^). The local gap enhancement 
factors A7A, are shown in Figure 2 (a) for various A, 
showing a modest increase of around 70% for A. = 1. The 
enhancement factor increases slightly with decreasing A 
but is essentially unchanged by modifications to phonon 
frequency and temperature for the parameter values used 
here. The temperature dependence of the gap was also cal- 
culated in Figure 2 (b). For very large temperatures, where 
ksT approaches A, there is a drop in the gap size. Below 
approximately 8,000 K (0.3^), this levels off, and the gap 
becomes relatively constant. 
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Figure 2 Modification of the BN bandgap. (a) The gap 

enhancement depends mainly on A., is weal<ly dependent on A and 
shows almost no change with Q. Calculations are made for A = t 
corresponding to a BN gap of 2 A = 4.66 eV, A = 1 .20t 
(2A = 5.6 eV), and A = 0.84t (2A = 3.92 eV, the tight binding fit 
from reference [1 0]). t=2.33 eV, hQ = 0.0]t = 23 meV, 
HQ = 0.03t = 70 meV, HQ = 0.06t = 140 meV, 
HQ = 0.09t = 210 meV, covering the full range of phonon 
frequencies in reference [1 3]. /cgT" = 0.01 1 (7=268K) and A. < 1 . (b) 
Variation of the gap with temperature, HQ = 0.05t =117 meV and 
k = 0.2. There is a weak temperature dependence due to the large A, 
consistent with the measurements in reference [1 5], with gap starting 
to close only for extremely high temperatures T > 8, 000 K, 
presumably above the melting point of the material. Red circles show 
the size of the gap, red lines are a guide to the eye. 



Conclusions 

A theory for the modification of BN band-gaps by interac- 
tion with phonons was presented here. It is of interest to 
make a comparison between the bandgaps of bulk h-BN, 
nanotubes, monolayer h-BN, and the theory presented 
here. Measured bandgaps of bulk h-BN are of between 
5.8 eV [15] and 5.971 eV [11], indicating that interaction 
between layers increases the bandgap, consistent with the 
theory here. The bulk gap is also higher than that for nan- 
otubes (5 eV) [16]. On the other hand. Song et al. [12] 
claim that the gap is reduced as BN thickness increases. 
The above discussion is presented with a warning that the 
theory requires that hopping between the substrate and 
the BN monolayer is small. Interlayer hopping will affect 
the bandwidth and bandgap, and the direct Coulomb 
interaction with strongly ionic substrates could also affect 
the band structure if the charge density at the surface of 
the substrate varies dramatically. 

It is also of interest to estimate the magnitude of the 
bandgap modification due to electron-phonon interac- 
tion in isolated monolayers of BN. Ab initio calculations 
have attempted to quantify the magnitude of the inter- 
action between electrons and acoustic phonons for small 
momentum excitations [17]. Extrapolating the interaction 
and taking a mean-field average (assuming mean momen- 
tum magnitude of ^7T/9a)y the electron-phonon coupling 
can be estimated as A. = {^n /9)^E\/2a^tMQ? , taking 
El = 3.66 eV from reference [17], M ^ 12.5 amu, a = 
2.5A. The mean energy of longitudinal acoustic phonons 
lies in the range of 50 to 75 meV, giving a range of A. = 0.05 
to 0.12, so the contribution of phonons to the bandgap is 
estimated as 3% to 7%. I would expect BN to have stronger 
interaction with optical phonons, since the pattern of dis- 
tortions around an electronic defect is consistent with 
optical modes (see Figure 1). 

The BN gap is too wide for digital applications. Recently, 
it has become possible to manufacture silicene, an atom- 
ically thick layer of silicon with similar properties to 
graphene [18], so it may be possible to make GaAs or AlP 
analogues to BN. Smaller gaps could be available from 
those materials, which might be used to create tunable 
bandgaps for atomically thick transistors. 
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